Projected free energies for polydisperse phase equilibria 
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A 'polydisperse' system has an infinite number of conserved densities. We give a rational proce- 
dure for projecting its infinite-dimensional free energy surface onto a subspace comprising a finite 
number of linear combinations of densities ('moments'), in which the phase behavior is then found 
as usual. If the excess free energy of the system depends only on the moments used, exact cloud, 
shadow and spinodal curves result; two- and multi-phase regions are approximate, but refinable 
indefinitely by adding extra moments. The approach is computationally robust and gives new geo- 
metrical insights into the thermodynamics of polydispersity. 
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The thermodynamics of mixtures of several chemical 
species is, since Gibbs, a well established subject (see 
e.g. [jjj). But many systems arising in nature and in in- 
dustry contain, for practical purposes, an infinite number 
of distinct, though similar, chemical species. Often these 
can be classified by a parameter, a, say, which could be 
the chain length in a polymeric system, or the particle 
size in a colloid; both are routinely treated as continuous 
variables. In other cases (see e.g. a is instead a 

parameter distinguishing species of continuously varying 
chemical properties. The thermodynamics of polydisper- 
sity (thus defined) is therefore of crucial interest to wide 
areas of science and technology. 

Standard thermodynamic procedures Q for construct- 
ing phase equilibria in a system of volume V containing 
M different species can be understood geometrically in 
terms of a free energy surface f(pj) (with / = F/V) in 
the M-dimensional space of density variables pj . Tangent 
planes to / define regions of coexistence, within which the 
free energy of the system is lowered by phase separation. 
The volumes of coexisting phases follow from the well- 
known 'lever rule' p|. Here 'surface' and 'plane' are used 
loosely, to denote manifolds of appropriate dimension. 
This procedure becomes unmanageable, both conceptu- 
ally and numerically, in the limit (M — > oo) of a polydis- 
perse system. There is now a separate conserved density 
p(a) for each value of a; the overall density of particles 
is p = J da p(a). The free energy surface is / = f[p(&)] 
which resides in an infinite dimensional space. Gibbs' 
rule allows the coexistence of arbitrarily many thermo- 
dynamic phases. 

Experimentally, one often restricts attention to the 
cloud- and shadow-curves (also referred to as dew/bubble 
curves). For a fixed 'shape' of polydispersity p(cr) = 
p(a)/p, these define as a function of the overall parti- 
cle density p (and temperature T) the onset of two-phase 
equilibrium (cloud-curve) and the density of the corre- 
sponding minority phase (its shadow); see e.g. 
Theoretical work has likewise focused on attempts to 
bring the problem into a more manageable form by some- 
how reducing its dimensionality P p|,p|-p7[ . In this Letter 
we propose a new and general method, whereby we con- 



struct from f[p(cr)] an (optimally) projected free energy 
surface in a reduced subspace of density variables. For 
these variables, we choose linear combinations of den- 
sities, the 'generalized moments' m, = J da Wi(a) p(a) of 
p{a), defined by certain weight functions Wi(a); these are 
ordinary (non- normalized) moments if Wi (a) = a 1 . 

The simplest imaginable case is where the free energy 
/ depends only on a finite set of K such moments: 



/ = f(m), 
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In coexisting phases one demands equality of particle 
chemical potentials, defined as p(a) = Sf/Sp(a) — 

J2i(df ' / dmi)wi(<j) = X«/^ w i( cr )i f° r au a - But tnis 
implies that all 'moment' chemical potentials, /ij = 
df/drrii, are likewise equal among phases. The osmotic 
pressures LT of all phases also must be equal; simple 
algebra establishes that —II = / — [lirrii which in- 
volves only the moments m,; and their chemical poten- 
tials pi. Finally, if the overall cr-distribution is p(°\a), 
and there are p coexisting phases with er-distributions 
p( a \cr), each occupying a fraction cf>^ of the total vol- 
ume (a = l...p), then conservation of particles im- 
plies the usual 'lever rule' (or material balance) among 
species: J2 a <t> ia) P (a) { a ) = P^V), Vcr. Multiplying this 
by a weight function Wi(a) and integrating over a shows 
that the lever rule also holds for the moments: 
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These results express the fact that any linear combination 
of conserved densities (a generalized moment) is itself a 
conserved density in thermodynamics. Therefore, if the 
free energy of the system depends only on K moments 
mi . . . m,K we can view these as the densities of K 'quasi- 
species' of particles, and construct the phase diagram 
via the usual construction of tangencies and the lever 
rule. Formally this has reduced the problem to finite di- 
mensionality by a projection, although this is trivial here 
because /, by construction, has no dependence on any 
variables other than the mi (i = 1 . . . K). 
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Of course, it is uncommon for the free energy / to 
obey (0). In particular, the 'ideal gas' (or, for poly- 
mers, Flory-Huggins) entropy term, in mixtures of many 
species, is definitely not of this form. On the other 
hand, in very many thermodynamic (especially mean 
field) models the free energy takes the form (fee = 1) 

/ = f{rm) + T J da p(a) [In (p(a)/R(a)) - 1] , (3) 

in which the excess free energy / does depend only 
on K moments. Examples include polydisperse hard 
spheres ||, polydisperse homo- and copolymers 
and van der Waals fluids with factorized interaction pa- 
rameters Note that, in the ideal gas term of ||), 
we have included a dimensional factor R(a) inside the 
logarithm: since the resulting contribution is linear in 
densities, this has no effect in rigorous thermodynamics. 
However, it will play a central role in our approach. 

In principle, the phase equilibria stemming from (|^) 
can be computed exactly by a finite algorithm. Specif- 
ically, the spinodal stability criterion involves a K- 
dimensional square matrix |ll]-[l4|] whereas calculation 
of p-phase equilibrium involves solution of (p— 1)(A'+ 1) 
strongly coupled nonlinear equations. This method has 
certainly proved useful p| -^,p|, ^0|Jl4| , but is cumbersome, 
particularly if one is interested mainly in cloud- and 
shadow-curves, rather than coexisting compositions deep 
within multi-phase regions ||M . Various ways of sim- 
plifying the procedure exist [p plf Jl5| -p^[| , but there has 
been, up to now, no systematic alternative to the full 
computation. Note also that the nonlinear phase equi- 
librium equations permit no simple geometrical interpre- 
tation or qualitative insight akin to the familar rules for 
constructing phase diagrams from the free energy surface 
of a finite mixture. 

Our method instead proceeds by deriving from (j3|) a 
'projected' free energy that depends only on a finite set 
of moments. We argue that the most important mo- 
ments to treat correctly are those that actually appear 
in the excess free energy f(rrii). Accordingly we divide 
the infinite-dimensional space of er-distributions into two 
orthogonal subspaces: a 'moment subspace', which con- 
tains all the degrees of freedom of p(a) that contribute to 
the moments m, (this subspace is spanned by the weight 
functions Wi(a)), and a 'transverse subspace' which con- 
tains all remaining degrees of freedom (as can be varied 
without affecting the chosen moments m^). Physically, 
it is reasonable to expect that these 'leftover' degrees of 
freedom play a relatively minor role in the phase equi- 
libria of the system, a view justified a posteriori below. 
Accordingly, we now allow violations of the lever rule, so 
long as these occur solely in the transverse space. The 
'transverse' degrees of freedom, instead of obeying the 
strict particle conservation laws, are chosen so as to min- 
imize the free energy: they are treated as 'annealed'. If, 



as assumed above, / = /(rrij) only depends on the mo- 
ments retained, this amounts to maximizing the entropy 
in (||), while holding fixed the values of the moments m,-. 

At this point, the factor R(a) in (0), which is immate- 
rial if all conservation laws are strictly obeyed, becomes 
central. Indeed, maximizing the entropy over all distri- 
butions p(a) at fixed moments to^ yields 

p{a) = R(a) exp (^T A * w * ( 4 ) 

where the Lagrange multipliers A^ are chosen to give 

mi — J da Wi(a) R(a) exp ^2 XtWiia)^ (5) 

The corresponding minimum value of / then defines our 
projected (i.e., annealed) free energy 



/pr(mi) = f + T 



to 



(G) 



In the last term, too = J da p(a) is the 'zeroth moment' 
which is identical to the overall particle density p de- 
fined previously. If this is among the moments used for 
the projection, the resulting linear term can be dropped; 
otherwise it must be retained (with Too now expressed as 
a function of the mi, via the A;). 

Our maximum entropy method yields a free energy 
/pr(wi) which only depends on the chosen set of mo- 
ments: i.e., ([|) is of the form (|l|) |l8|. A finite dimen- 
sional phase diagram can now be constructed from it ac- 
cording to the usual rules. Obviously, though, the results 
now depend on R(a) which is formally a 'prior distribu- 
tion' for the entropy maximization. To understand its 
thermodynamic role, we recall that our projected free en- 
ergy /pr(wi) was constructed as the minimum of f[p{a)] 
at fixed to,; that is, / pr is the lower envelope of the pro- 
jection of / onto the moment subspace. Crucially, the 
shape of this envelope depends on how, by choosing a 
particular prior distribution R(a), we 'tilt' the infinite- 
dimensional free energy surface before projecting it. 

To find the optimum choice of prior, we note that 
R(a) serves physically to determine which distributions 
p(a) lie within the maximum-entropy family (|j) that 
the annealed system can have. Typically, one is in- 
terested in a system where a fixed overall 'parent' (or 
'feed') distribution p^> (a) becomes subject to separation 
into various phases. In such circumstances, we should 
generally choose this parent distribution as our prior, 
R(a) = p(°'(a), thereby guaranteeing that it is contained 
within the family (||). Having done this, we note that the 
annealing procedure will be exactly valid, to whatever 
extent the cr-distributions actually arising in the various 
coexisting phases of the system under study are mem- 
bers of the family (^). (This statement of exactness, and 
similar ones below, of course hold only if (0) is valid.) 
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In fact, the condition just described does hold when- 
ever all but one of a set of coexisting phases are of in- 
finitesimal volume compared to the majority phase. This 
is because the er-distribution, p(°'(cr), of the majority 
phase is negligibly perturbed, whereas that in each mi- 
nority phase differs from this by an exponential Gibbs- 
Boltzmann factor, of exactly the form required for (|4|). 
Accordingly, our projection method yields exact cloud- 
curves and shadow-curves. By the same argument, crit- 
ical points (which in fact lie at the intersection of these 
two curves) are exactly determined. Moreover, all spin- 
odals are also found exactly by our annealing method. 
For, at a spinodal, there exists an instability direction 
(in the full space) along which the curvature of the free 
energy vanishes; in all other directions / has positive 
curvature. One can show that such an instability direc- 
tion always connects neighboring distributions within the 
same maximum entropy family (Q), and hence that only 
the free energy of such distributions (i.e., the projected 
free energy with the parental R(a)) is needed to calculate 
spinodals. The geometrical interpretation of this result, 
and also proofs of it and the others stated above, will be 
given elsewhere Q. 

The method does, however, give only approximate re- 
sults for coexistences involving finite amounts of differ- 
ent phases. This is because linear combinations of dif- 
ferent cr-distributions obeying ([!]), corresponding to two 
(or more) phases arising from the same parent (p^ (u) — 
-R(cr)) do not necessarily add to recover the parent distri- 
bution itself. Moreover, according to Gibbs' phase rule, 
a projected free energy depending on n moments will 
not normally predict more than n + 1 coexisting phases, 
whereas a polydisperse system can in principle separate 
into an arbitrary number of phases. Both of these short- 
comings can be overcome by systematically including ad- 
ditional moments within the annealing procedure. (The 
above exact results are unaffected, because these do not 
exclude a null dependence of / on certain of the m,.) In- 
deed, by adding further moments one can indefinitely ex- 
pand the maximum-entropy family (^) of cr-distributions, 
thereby approaching with increasing precision the actual 
distributions in all phases present; this yields phase di- 
agrams of ever-refined accuracy. How quickly conver- 
gence to the exact results occurs depends on the choice 
of weights functions for the additional moments; this will 
be quantified elsewhere |HJ . 

To demonstrate the power of our approach, we consider 
a specific example. This is a simplified model of chemical 
fractionation, in which one considers species of continu- 
ously variable chemical character (such as aromaticity) 
governed by a parameter a between and 1. We sup- 
pose that the interaction energy between species varies 
as (<t — er') 2 , so that the most different species repel each 
other most strongly. For simplicity we take a molten sys- 
tem, choosing volume units so that the overall density is 
constrained as fdap(a) — 1. Within a mean-field treat - 
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FIG. 1. Coexistence curves for a parent distribution with 
mf = 0.6. Shown are the values of mi of the coexisting 
phases; horizontal lines guide the eye where new phases ap- 
pear. Curves are labeled by n, the number of moments re- 
tained in the projected free energy. Predictions for n = 10 
are indistinguishable from an exact calculation (in bold). 

ment, the system is then described by a free energy of 
the form (JsJ) , with an excess free energy (in units of k^T) 
of / = — x m ii tui = Jdaap(a) (up to irrelevant terms 
linear in p{a))- This model differs only by a reseating of 
parameters (with powers of polymer molecular weight) 
from the Flory-Huggins treatment of random AB copoly- 
mers, in which \ is the usual interaction parameter and 
a is the proportion of A monomers in a chain ]^-^| . The 
model should show fractionation into an ever-increasing 
number of phases as \ is increased. It is therefore an in- 
teresting test case for our projection approach (and the 
method of adding further moments), yet simple enough 
for exact phase equilibrium calculations to remain feasi- 
ble, allowing detailed comparisons to be made. 

We consider phase separation from parent phases with 
a distributions of the form p(a) oc exp(Acr) (for < a < 
1); A is thereby fixed in terms of the parental mi = . 
Fig [l] shows the exact coexistence curve for = 0.6, 
along with the predictions from our projected free energy 
with n moments (m., = J da a 1 p{a) , i = 1 . . .n) retained. 

Comparable results are found for other m[ . Even for 
the minimal set of moments (n = 1) the point where 
phase separation first occurs on increasing \ is predicted 
correctly (this is a cloud point for the given parent). 
As more moments are added, the annealed coexistence 
curves approach the exact one to higher and higher pre- 
cision ^0|. As expected, the precision decreases at high 
X, where fractionated phases proliferate; in this region, 
the number of coexisting phases predicted by the projec- 
tion method increases with n. However, it is not always 
equal to n + 1 , as one might expect from a naive use of 
Gibbs' phase rule; three-phase coexistence, for example, 
is first predicted for n = 4 pl[|. Note that the stability 
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of the results to the addition of extra moments provides, 
in this example, a good test of convergence on the coex- 
istence curves. 

For the computational implementation of both the 
annealing and the exact method we used a Newton- 
Raphson nonlinear equation solver. The annealed cal- 
culation turned out to be significantly more robust with 
respect to the choice of initial values, the size of x m - 
crements etc. due to an effective decoupling of the equa- 
tions: Equality of chemical potentials is achieved using 
the moments contained in the excess free energy, while 
the lever rule is satisfied (increasingly accurately) using 
the remaining moments |19f |. This advantage should be 
much more pronounced in more complex cases, as should 
savings in computer time (which are modest in our simple 
example) . 

With exact results for cloud- and shadow-curves, crit- 
ical points and spinodals, as well as refinably accurate 
coexistence curves and multi-phase regions, our anneal- 
ing method allows rapid and accurate computation of the 
phase behavior of many polydisperse systems. Moreover, 
by establishing the link to a projected free energy f w {m,i) 
as a function of a finite set of conserved densities m^, it 
restores to the problem much of the geometrical interpre- 
tation and insight (as well as the computational method- 
ology) associated with phase diagrams for finite mixtures. 
This contrasts with procedures commonly used for sys- 
tems in which the excess free energy involves a finite set 
of moments <j3j> p|-p|,p|pO[ . Some previous approxima- 
tions to that problem have used (generalized) moments 
as coordinates; see e.g. (Tl|,|l|,^|Jn|j2|]. Our annealing 
method provides a rational basis for these methods and, 
by a careful choice of prior, guarantees that many prop- 
erties of interest are found exactly. 

Finally, our method may extend to models for which 
the excess free energy cannot be written directly in terms 
of a finite number of moments as in (J3J) - For exam- 
ple, many mean-field theories correspond to a variational 
minimization of the free energy: F < (E}q — TSq, where 
subscript refers to a trial Hamiltonian In such a 
case, one might choose to first make a physically moti- 
vated decision about which (and how many) moments rrii 
to keep, and then include among the variational param- 
eters the annealed "transverse" degrees of freedom. This 
would lead directly to a mean-field estimate of the pro- 
jected free energy without assuming Eq. (||). Note that 
a good choice of prior R(o~) will again be important. Al- 
though no exact results can be guaranteed, this approach 
may form a promising basis for future developments. 
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